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Abstract 

We investigate the influence of the substrate on the vibrational properties of graphene, 
comparing graphene on hexagonal boron nitride (hBN) with graphene on SiC>2 by spatially 
resolved confocal Raman spectroscopy. By studying the G-line we show that the average 
doping level and local doping domain fluctuations are significantly suppressed in graphene on 
hBN with respect to graphene on S1O2. However, in contrast to the G-line, the 2D-line of 
graphene on hBN shifts up in frequency compared to the one of graphene on Si02- We show 
*To whom correspondence should be addressed 
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that this effect is due to a reduction of the Kohn anomaly at K through an enhanced screening 
by the dielectric substrate. We prove our theory to be consistent with Raman measurements on 
graphene surrounded by hBN (stronger screening) and recent findings on suspended graphene 
(no external screening). 

Graphene, a monoatomic carbon membrane with unique electronic properties^ is a promis- 
ing candidate for flexible electronics, high frequency applications and spintronics. 3 However, 
graphene's ultimate surface-to-volume ratio makes the environment, in particular the substrate 
material, have a pronounced influence onto its intrinsic properties. For example, S1O2, the most 
common substrate material, exhibits surface roughness, dangling bonds and charge traps which in- 
troduce ripples, disorder, 4 and doping domain fluctuations. 5 This limits carrier mobilities and the 
operation of graphene devices .^Therefore alternative substrates are required to overcome these 
limitations. Hexagonal boron nitride (hBN) has been identified as a very promising candidate 
A large (indirect) band gap, a lattice mismatch to graphite of less than 2%, and the absence of dan- 
gling bonds makes this atomically flat material a valuable and promising insulating counter part to 
grapheneJ 12 * 13 ! It has been shown that graphene can be successfully transferred to ultra-thin hBN 
flakes leading to improved electronic transport properties compared to graphene on S1O2. lllMore- 
over, scanning tunneling microscopy experiments have shown that the sizes of individual electron- 
hole puddles are significantly increased while the disorder potential is reduced by roughly a factor 
10.™ 1 Over the last years Raman spectroscopy has proven to be a powerful tool for characterizing 
graphene and studying its physical properties. For example, this technique has been successfully 
used (i) to distinguish single-layer graphene from few-layer graphene and graphiteJ^E^ (ii) to 
monitor doping levels,'^' (iii) to study short range disorder and edge properties^ and (iv) to in- 
vestigate suspende and nanostructured graphene. I^Here we present spatially resolved confocal 
Raman spectroscopy measurements of graphene on hexagonal boron nitride substrates, which are 
compared with measurements of graphene on S1O2. We show that the average doping level and 
local doping domain fluctuations are significantly suppressed in graphene on hBN. This is visible 
in a redshift of the G-line which has also been observed for suspended graphene.^ The 2D-line of 
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Figure 1: (color online) (a) Schematic illustration of a graphene sample with incident laser, (b-c) 
Optical microscope images of graphene flakes partly resting on hBN (blue) and Si02- The flake 
consists of different regions of single-layer, bilayer and few-layer graphene. (d) Raman spectrum 
of single-layer graphene resting on hBN. Inset: Region around the hBN peak at sites marked in 
panel c. (e-f) G (e) and 2D (f) peak of single-layer graphene on hBN (marker A) and Si02 (marker 
D). 



graphene on hBN, however, displays an opposite behavior to that of suspended graphene. For sus- 
pended graphene, the 2D-line is red-shifted (compared to graphene on SiO^)^ while we observe a 
significant blue-shift with respect to graphene on SiC>2. We show in this letter that this difference 
is due to the presence/absence of dielectric screening which has an influence on the frequencies of 
the highest optical phonon branch. 

Graphene and hBN flakes are prepared by micromechanical cleavage. While hBN is directly 
deposited on a Si02/Si ++ substrate, graphene is prepared on top of a polymer stack consisting of 
a water-soluble polymer [100 nm Polyvinylalcohol (PVA)] and a water resistant polymer [270 nm 
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Polymethyl methacrylate (PMMA)] allowing the transfer process described in detail in Ref. 8. Be- 
fore depositing the graphitic flakes on top of the hBN, Raman spectroscopy has been used to 
identify and select individual single-layer graphene flakes J^E^ A schematic illustration of a sam- 
ple is shown in Figure la and optical microscope images of some fabricated samples are shown in 
Figures lb-c. 

The Raman data are recorded by using a laser excitation of 532 nm (£^=2.33 eV) through 
a single-mode optical fiber whose spot size is limited by diffraction. A long working distance 
focusing lens with numerical aperture of 0.85 is used to obtain a spot size of approximately 400 nm. 
We used a laser power below 1 mW such that heating effects can be neglected.^ 

A typical Raman spectrum of a graphene flake resting on hBN is shown in Figure Id. The 
Raman spectrum shows the prominent G-line around 1582 cm -1 as well as the single Lorentzian 
shaped 2D-line around 2675 cm -1 as expected for graphene. A third prominent and sharp peak 
arises around 1365 cm -1 , which can be attributed to the Raman active LO-phonon in hBN.^It is 
important to distinguish this peak from the defect induced D-line potentially appearing at around 
1345 cm^ 1 .^ Therefore, Raman spectra have been acquired at edge regions of the graphene flake 
where defects are known to appear. 17 The insets in [TJl show corresponding Raman spectra at 
different positions marked and labeled infTb. The data recorded at the edge (B) shows a second peak 
arising at around 1345 cm -1 which is not visible in the bulk region and can be clearly distinguished 
from the one at 1365 cm -1 . As shown in inset C of [TJi, the D-line also appears in regions of 
substrate transition, i.e. at the edge of the underlying hBN flake. This can be attributed to local 
bending of the graphene flake induced by the level difference of the two substrates. The sp 2+ri 
hybrid orbitals necessary to bend the graphene layer cause short range scattering leading to a D 
peak in the Raman spectrum. The D line is not visible in regions of graphene away from the 
edges, neither on hBN nor on Si02- We conclude that the transfer technique used in the fabrication 
process does not induce a significant amount of defects in the graphene lattice. 

In order to show the substrate dependence of the Raman lines of graphene, we compare in 
Figures le-f the G and 2D-lines of a flake that partially rests on hBN and partially on Si02- The G 
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Figure 2: (color online) (a) 2D FWHM Raman image of the sample shown in [lp. Where no 2D 
peak was found the hBN peak FWHM was plotted making the hBN flake visible (highlighted by 
white dotted lines), (b) Raman map of the 2D FWHM of the sample shown in Figure lb (see white 
dashed box in Figure lb). Dotted line marks the hBN substrate, (c) Raman map of the same sample 
but integrated intensity of the hBN peak, (d) Line cut along the solid arrow (path P) in panel a. 

peak of graphene on hBN and the one on Si02 differ significantly in their position: the G peak on 
SiC>2 is centered at 1586.5 cm -1 while the one on hBN is centered at 1582.8 cm -1 . This downshift 
can be attributed to reduced doping, which also is consistent with the increase of the full width at 
half maximum (FWHM) of the G peak of graphene on hBN.^The FWHM is 12.2 cm 1 on Si0 2 
and 16.7 cm -1 on hBN. Also the 2D peak shows a substrate dependence of the position which is 
2674.0 cm -1 on SiC>2 and 2681.6 cm -1 on hBN. The substrate dependence of the FWHM is even 
more significant, being 36.4 cm -1 on Si02 and 28.1 cm -1 on hBN. 

The peak-shifts and changes in the FWHM can not only be seen in individual Raman spectra, 
but also appear spatially resolved in two dimensional Raman maps. A Raman map of the 2D 
FWHM of the sample presented in [TJ; is shown in [2^. One can identify three single layer regions 
with a FWHM below 40 cm -1 , two resting on hBN and one resting on SiC>2 with a small region 
also resting on hBN. A line cut in this substrate transition region shown in [2]i reveals a locally 
resolved difference of the FWHM of around 8 cm -1 . [2p- shows the 2D FWHM map of the sample 
previously shown in[T[) (left panel). There is also a substrate dependency visible. A Raman map 



of the integrated peak intensity between 1360 and 1370 cm is shown in[2p. The bright area has 
exactly the same shape as the underlying hBN flake in the optical picture and verifies the attribution 
of this peak to the hBN mode. 
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Figure 3: (color online) Statistical evaluation of single-layer regions of the Raman map shown in 
[2^ in terms of G peak position (a), 2D peak position (b), intensity ratio between G and 2D peak 
(c), and FWHM of the 2D peak (d). 



To relate the quantitative results of the measurements to physical properties like charge carrier 
density fluctuations and to dispose of statistical fluctuations, it is necessary to evaluate a large num- 
ber of spectra statistically. The statistical distribution of the G peak position of the flake of [TJ;, only 
considering single-layer regions, is plotted in [3^. These data show a clear distinction between the 
different substrates. To obtain estimates of the statistical parameters the distribution is considered 
to be approximately Gaussian. The mean value of the distribution on hBN, /ihBN = 1583.1 cm -1 , 
is red shifted with respect to the one on SiC>2, jUsi0 2 =i 585.9 cm -1 , a notable deviation by almost 
three wavenumbers. Even more pronounced is the difference in the standard deviation which is a 
measure for the G peak fluctuations. Being OhBN=0.7 cm -1 on hBN, it is almost four times smaller 
than on Si02 where 0siO 2 =2-7 cm -1 . 

Previous Raman measurements with gated graphene flakes on SiC>2 have demonstrated a de- 
pendency of the G peak position on the charge carrier density J^^S A non-adiabatic theory was 
established to calculate the G peak shift in terms of charge carrier densities.! 19 ! 27 ^ By using a fi- 
nite temperature of 295 K and an intrinsic G peak position of 1582.5 cm -1 according toj^we 
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obtain a charge carrier density of 1.8 xlO 12 cm" 2 on Si02 and 9xlO n cm -2 on hBN, meaning 
the overall doping of the investigated single-layer flake on hBN is reduced by a factor of two with 
respect to a region of the very same flake resting on Si02- Furthermore, using the standard devi- 
ation of the G peak shift distribution to quantify the charge fluctuations, one obtains a fluctuation 
of 1.6x 10 12 cm -2 on SiC>2 and a fluctuation of 6x 10 11 cm -2 on hBN. This difference by almost 
a factor of three indicates a significant reduction in doping domain fluctuations and hence a re- 
duction of the disorder potential in the single-layer graphene. Comparing these results with the 
data obtained by scanning tunneling spectroscopy } w ^ one will notice a difference by one to two 
orders of magnitude. However, the experiments using scanning tunneling spectroscopy, besides 
being carried out in a low temperature controlled environment, acquire data with an atomic res- 
olution from an area of 100 nm 2 which is on the scale of individual charge puddles, the Raman 
setup with a spatial resolution of around 500 nm is capable of measuring on a micrometer scale 
averaging over a large number of charge puddles. For instance, the areas used to acquire the data 
on this sample are around 50 /im 2 on Si02 and around 100 /im 2 on hBN. 

As shown in previous works, the reduced doping fluctuations are also manifested in a reduced 
ratio between the integrated peak intensities of the G and 2D peak.^EH j^g statistical distribution 
depicted in [3]: shows a clear substrate dependency of this ratio and is in qualitative agreement with 
the G peak evaluation. 

Experimental data of the statistics of the FWHM of the 2D-line is provided in [3ji. The mean 
values for the two substrates are clearly different and also the standard deviations differ. Con- 
sistent with the two dimensional Raman map shown before, the 2D peak width of the regions 
on hBN (/ij 1 BN = 25.2 cm -1 ) is significantly smaller by almost 10 cm -1 than the one on Si02 
(/isio 2 =34.5 cm -1 ). Also the fluctuations of the hBN data, ChBN=l-3 cm -1 , are significantly 
smaller than the ones of the SiC>2 data, 0sio 2 = 2.7 cm -1 . So far, the FWHM of the 2D peak 
was not considered to be doping dependent. Measurements of gated graphene on Si02 showed 
no significant dependence on the charge carrier density.^ Hence, a reduced doping alone cannot 
explain this substrate dependence. A similar reduced line- width of the 2D-line (23 cm -1 ) was 
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observed for suspended graphene. 21 We assume that the increased FWHM for graphite on Si02 
is due to the substrate roughness and the presence of impurities which gives rise to an enhanced 
phonon scattering and thus a smaller life-time of the phonons that are excited during the 2D-Raman 
process. 

We now turn to the discussion of the 2D-line position whose statistical distribution on hBN and 
on SiC>2 is shown in [3}). While the position of the G-line can be directly related to the presence 
or absence of residual charging due to impurities on or in the substrate, the interpretation of the 
2D-line shift is more subtle. Small charge densities (< 4xl0 12 cm~ 2 ) lead to shifts of the 2D 
line by at most 2 cm -1 (Ref. 5). Furthermore, the 2D-line positions of suspended graphene^ and 
graphene on hBN differ by almost 10 cm -1 even though both systems are mostly free of charging 
impurities, as demonstrated by the coincidence of the G-line positions. Thus, we discard charging 
as the origin for the 2D-line shift. 

Our interpretation of the 2D-line shift is based on the double-resonance Raman model of Thom- 
sen and Reich. 29 The model successfully describes the D and 2D dispersion as a function of laser 
energy as well as the splitting of the 2D-line for bilayer, and few-layer graphene,! 16 * 17 ! 30 ! provided 
that renormalization of the highest optical phonon branch (HOB) due to electron-correlation effects 
is properly taken into account. 31 According to the double-resonance model, the 2D-line dispersion 
is proportional to the slope of the HOB between the high-symmetry points K and M, and inversely 
proportional to the slope of the % bands around K. In a first step, we thus calculated the electronic 
band structure and the phonon-dispersion of graphene on hBN using standard density-functional 
theory in the local density approximation (DFT-LDA).^ We chose the most stable configuration 
where one carbon atom is on top of a boron atom and the other carbon atom in the hollow site 
of hBN.^1 In the electronic band-structure a small gap of 53 meV opens at K but further away 
from K, the % -bands remain almost unchanged.^ In the phonon dispersion, calculated by density- 
functional perturbation theory (DFPT), 35 there is only a small change of the HOB in the immediate 
neighborhood around K. There, one observes a slight smearing of the Kohn anomaly (due to the 
very small band-gap opening), manifest in an up-shift of the HOB by about 3 cm -1 . Everywhere 
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Figure 4: (color online) (a) Geometry employed for the calculation of the electron-phonon cou- 
pling in graphene surrounded by hBN. (b) Schematic figure of the differences in the highest-optical 
phonon branch around K between suspended graphene and graphene surrounded by hBN. (c) Ra- 
man spectrum of graphene surrounded on both sides by multi-layer hBN. The sample is shown in 
the left inset. The darker-blue area is the region of graphene on an underlying hBN flake. In the 
light-blue area, an additional hBN layer is deposited on top of the graphene layers. The spot where 
the spectrum was measured is marked by the black cross and the scale bar is 5 /im. (d) Raman 
2D-line of graphene on SiC>2, graphene on hBN, and graphene surrounded by hBN. 



else between K and M, in particular in the wave-vector range that is sampled in Raman exper- 
iments, the upshift of the HOB is less than 1 cm -1 . We conclude that the pure "mechanical" 
interaction alone between graphene and the hBN substrate cannot explain the blue-shift of the 
Raman ID-line. 1 ^ 1 

In recent work, it was shown by calculations on the level of the GW- approximation 38 that elec- 
tronic correlation beyond DFT-LDA influences both the slope of the 7T-bands of graphene 39 l 4 ° l and 
the slope of the highest-optical phonon branch around K. 31 The electron-electron interaction de- 
pends on the electronic screening by the environment. Therefore, we expect that correlation effects 
in graphene will be reduced by a dielectric substrate. Although SiC>2 and hBN have both roughly 
the same dielectric constant the coupling between graphene and ultra-flat hBN is significantly in- 
creased compared to graphene on rough SiC>2 (see also illustration in Figure la). This different 
dielectric environment will have consequences for the Fermi velocity and for the electron-phonon 
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Table 1: Calculated band-gap opening (for a displacement d = 0.0053 A) and electron-phonon 
coupling of the highest optical (A[) phonon at K. Comparison of LDA and GW calculations. 
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coupling. In order to verify this hypothesis, we have performed GW-calculations on isolated (sus- 
pended) graphene and on graphene surrounded by two layers of hBN. The periodic geometry that 
we used in our calculations is shown in [4^. For simplicity, we have symmetrized the unit cell by 
choosing an ABC stacking sequence for the three layers. This ensures that the two carbon atoms 
in the unit-cell are equivalent, each with a boron atom on one side and a hollow site on the other 
side. Due to the symmetry, the linear crossing of the 7T -bands is preserved and we can use the same 
strategy as in Ref. 31 for the calculation of the electron-phonon coupling at the high-symmetry 
point K: In the y/3 x y/3 supercell, the atoms are displaced by a small distance (d = 0.01 atomic 
units) from their equilibrium position along the eigenvector of the HOB (see Figure 3b of Ref. 31). 
The squared electron-phonon (e-ph) coupling - which determines the slope of the HOB around 
- is then obtained as 

. \ 1 /A£ K ^ 2 



where AEj^ is the induced energy gap between the n* and % bands at K. 

In Table 1, we present the results of our calculations. 42 On the LDA level, the e-ph coupling 
is 3.7% weaker for graphene surrounded by hBN than for pure (suspended) graphene. This differ- 
ence is increased to 8% on the level of the GW approximation. Obviously, the increased screening 
by the BN substrate reduces the gap opening and thus the e-ph coupling. 44 Since the slope of the 
HOB around K is proportional to the e-ph coupling 41 and since the frequency of the HOB far 
away from K is almost independent of screening effects,^ a reduction of the e-ph coupling leads 
to an increase of the HOB frequency at and around K. This is demonstrated in [4J3. In order to 
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quantitatively understand this effect, it would be desirable to calculate the exact phonon disper- 
sion relation including correlation effects. Due to the complexity of total-energy GW calculations, 
this is currently not feasible. A workaround consists in the use of the hybrid B3LYP functional^ 
where we adjust the parameters in order to mimic the two different screening values for "pure" and 
"sandwiched" graphene.^ Using these functionals, we obtain that the HOB shifts (from "pure" 
to "sandwiched" graphene) by +16 cm -1 at K, by +2.4 cm -1 at M, and by +6.8 cm -1 half-way 
between K and M (which corresponds to a excitation energy of 2.8 eV). At the same time, the 
frequency at T remains almost unchanged (-0.4 cm -1 ). For the laser energy of 2.33 eV, linear in- 
terpolation yields a phonon shift of +8.3 cm -1 . The corresponding 2D-line shift (where 2 phonons 
are excited/absorbed) would be then +16.6 cm -1 . These calculations are not meant to provide 
absolute numbers for the 2D-line shift but they demonstrate that the experimentally measured 2D- 
line shift qualitatively agrees with the shift that is induced by the enhanced dielectric screening for 
graphene on hBN.^1 

So far, we have compared a calculation for graphene surrounded by hBN with experiments 
where graphene is deposited on top of a hBN flake. This motivated us to perform Raman mea- 
surements on graphene surrounded by hBN on both sides. We achieved this by the deposition of 
an additional multi-layer hBN flake on top of one of our graphene on hBN samples. The result- 
ing Raman spectrum is shown in [4]:. The spectrum displays three prominent peaks: the peak at 
1367 cm -1 is the E2g2" m °de of hBN. The G-line of graphene at 1583.7 cm -1 is approximately in 
the same position as the G-line of graphene with hBN on one side. This evidences that only a few 
additional charge impurities are added during the deposition of the top-hBN flake. In contrast to 
the G-line, the 2D-line of surrounded graphene around 2687.4 cm -1 is considerably blue-shifted 
by 7 cm -1 compared to the 2D-line of one-sided graphene on hBN as shown in Figure 4d. With 
respect to suspended graphene (2673.5 cm -1 ) 21 the blue-shift of the 2D-line of graphene on hBN 
roughly doubles when a second hBN-layer is added on top. This is another indication that the 
blue- shift has its origin in the screening dependence of the HOB between K and M. 

In conclusion, our Raman measurements confirm that hBN is a high-quality insulating sub- 
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strate for graphene with strongly reduced impurity charging. This is evidenced by the red-shift of 
the G peak with respect to graphene on a standard Si02 substrate. In contrast to the G-line, the 
2D-line of graphene on hBN is blue-shifted. We have shown that this change in the frequency 
of the highest-optical branch is not a consequence of a direct (mechanical) interaction between 
graphene and its substrate. It is rather an indirect, electronic, effect mediated by the influence of 
the dielectric screening on the electronic structure of graphene. Usually, in metallic systems, the 
phonon-frequencies close to a Kohn anomaly depend only on the internal screening of the (bulk) 
material. Layered graphene is an example where (close to the Kohn anomaly at K), the frequencies 
of the highest optical branch depend on the external screening by the local dielectric environment. 
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